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' Abstract 

S. 

This article deals with transport properties of one dimensional Brownian diffusion 
under the influence of a correlated quenched random force, distributed as a two-level 
Poisson process. We find in particular that large time scaling laws of the position 
of the Brownian particle are analogous to the uncorrelated case. We discuss also the 
probability distribution of the stationary flux going through a sample between two 
prescribed concentrations, which differs significantly from the uncorrelated case. 
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Disorder significantly influences diffusive transport phenomena [1] [2] [3]. Indeed, 
strong enough disorder does not induce a simple renormalization for transport coeffi- 
cients of the corresponding pure system, but usually generates "anomalous diffusion" 
effects. The case of spatially uncorrelated disordered media is now pretty well under- 
stood. The effects of spatial correlations on the dynamics have first been discussed for 
the directed random walk problem [4]. More recently, the response to an external ap- 
plied force has been studied for the one dimensional case of a Brownian particle diffusing 
under the influence of a quenched random force {F(x)} [5]. For a sample characterized 
by some particular realization of the stochastic process {F(x)}, the diffusion is defined 
by the following Fokker-Planck equation for the probability density P(x,t |xo,0) 

dP d ( dP \ 

In this article, we analyse in details transport properties of this model in the case where 
{F(x)} is distributed as a two- level Poisson process. 

Let us define more precisely the model of disorder we consider and some of its 
properties. We assume that the quenched random force F(x) takes alternatively a 
positive value 0o > and a negative value — 0i < on intervals whose lengths are in- 
dependent random variables distributed according to the following probability densities 
respectively 

f (l) = 9(1) n e~ n ° l 

2 

f x (h) = 9(h) m e~ nih 

The parameters ^ and ^- are respectively the mean length of intervals {F(x) = 0o} 
and the mean length of intervals {F(x) = — 0i}. This choice of exponential distributions 
for fo and j\ is in fact the only one that makes the process {F(x)} Markovian. This 
property enables us to write differential equations for the probability po(x) to have 
F(x) — 0o an d the probability p±(x) = 1 — po(x) to have F(x) = — 0i. 

dpo 




= -n po + n\p\ =m- (n + ni)p 
= ~n\p\ + n p = n - (n + n ± )p ± 



(3) 



dx 

The mean value F and the two point correlation function G(x) of the process {F(x)} 
read 

n\ no 

F =< >= 0o : 01 : 

n + ni n + ni 

( 4 ) 

G(X) =< 0(X)0(O) > - < > 2 = n0Ul (0o + 1 )2 e -(-o+n 1 )|x| 

(n J rn 1 ) 1 

The corresponding random potential U(x) = — J x F(y)dy seen by the Brownian 
particle presents an alternance of positive and negative slopes of random lengths as 
sketched on Figure 1. 
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The fundamental random variable associated with classical diffusion under the 
action of a quenched random force presenting a stricly positive mean < F(x) >= F > 0, 
is the exponential functional of the random potential U (x) 



/•OO /"OO 

t oq = dx eP U ^ x > = dxe Jo (5) 
Jo Jo 



-(3 [ F(y)dy 
Jo 



Indeed, the probability distribution of the functional determines the large time 
anomalous behaviour of the Brownian particle position [3]. In particular, the velocity 
defined for each sample as 



V = lim — / dx xP(x,t \xq,0) (6) 

t^oo dt J_ 00 

is a self-averaging quantity [6] (for any homogenous random potential presenting only 
short-range correlations) inversely proportional to the first moment of [3] 

V=^- (7) 

<r ao > 

When the quenched random force is distributed with the Gaussian measure 

-i- / dx [F(x) - F ] 2 
VF(x) e 2a J (8) 

the probability distribution 7 3 00 (r) of the functional reads in terms of a = and 
V = it > [3] 

a ( 1 -— 1 



Voo{t) = —— — e dT oc — — (9) 

The algebraic decay at large r explains all the dynamical phases transitions between 
different anomalous behaviours known for this model [3]. In particular, Eq (7) implies 
that the value fx = 1 separates a phase of vanishing velocity V = for < \i < 1, and 
a phase of finite velocity V > for \x > 1. Another interesting physical quantity is the 
stationnary current Jn which goes through a disordered sample of length iV between 
a fixed concentration Pq and a trap described by the boundary condition Pjy = [7] 
[8]. In the limit iV — > oo, the stationary flux Joo is simply a random variable inversely 
proportional to the functional 

Joo = ^ (10) 

Too 

Note that unlike the velocity V, this flux is not a self-averaging quantity, but must be 
described by its full probability distribution. 
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It is very convenient to introduce the more general functional 



[ d yeW(y)-V( X )\ = f dye ^L ^ („) 

J X J X 

and to consider the random variable as the limit of this process as x — > — oo 

lim r(x,b) (12) 



' oo 

x — ► — oo 



The evolution of the functional t(x, b) is governed by the the stochastic differential 
equation 

^-=f3F(x)r(x,b)-l (13) 

The stochastic term F(x) appears multiplicatively, so that the fluctuations of the ran- 
dom force are coupled to the values taken by the random process r(x, b) itself. When 
the random force {F} is distributed as a white noise (8), the multiplicative stochastic 
process t(x, b) can be in fact be related to Brownian motion on a surface of constant 
negative curvature [9]. 

We now compute the probability distribution of the functional when the 
quenched random force {F(x)} is a two-level Poisson process, and compare it with 
the result (9) for the white noise case. The random variable only exists if the force 
mean value F is strictly positive 

(no + ni)F = (porn — (pino > that is ^ — > (14) 

01 00 

that we assume from now on. Note that the random variable t(x, b) remains confined 
in an interval [r m i n (x, 6), r max (x, b)] depending on the length (b — x) 



, 1 _ p-l3<t>o(b-x) 

T min (x,b) = j b x dye-^y~^ = - 



Tmax(x,b) = gdy e+^ifo"*) 



/30o 

e /30i(6-x) _ 2 



(15) 



This is very different from the white noise case, where the random force is not bounded. 
Even in the limit where the interval length (b — x) tends to infinity, the support of the 
random variable is not [0, +oo], but [^^, +oo]. 

The two-level Poisson process {F} is Markovian; so according to the local evolution 
equation (13), the coupled process {F(x), t(x)} is still Markovian. Let us define the 
joint laws 

Pq(t, x)dr = Prob< F(x) = 0o and t(x, b) E [r, r + dr] 
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Pi(r,x)dr = Prob|F(x) = -0i and r(x,b) G [r, r + dr] } (16) 
They evolve according to the two coupled Master equations 



BP B 

' ° =- 7r [(l-f3MPo]-n Po + n 1 P 1 



Bx Br 

BP B 

' -— [(1 + /30it)Px] + n P - niPi 



(17) 



<9x <9r 



Stationary distributions Po(t) et Pi(r) in the limit x — > — oo are therefore the solutions 
of the system 

■f [((3<P t - 1)P ] - n P + n 1 P 1 = 



dr 

— [(pfar + l)Pj - n P + reiPi = 



(18) 



respectively normalized on the interval [^^, +oo] by 

1 P ( r )rfr = ^^ and /i P x { T )dr = — ^— (19) 



T • 11 

It is convenient to set To = — — , T\ = — — , v§ = n ro et v\ = n\T\, and 

p(Po P<Pi 

v=vi -* = n*-£) >0 (20) 

which is strictly positive according to the hypothesis Po > (14). The solutions of (18) 
(19) read, using Heaviside function 6 , 

P ( T ) = Ag(T-7a) L_i^ where A = ^ r^^g^ (r + n)^ 

P 1 (r) = P^-ro) ( ) T + " ri r ^; i where B = ^t^^ (t + n)— 

(21) 

Moments of order of these distributions diverge as soon as /c > z/ = (z/i — z/o)- They 
read otherwise in polynomial forms as 



(22) 



/ °°dTT*P (r) =^t~ r k fU % 1 + ^-^i;1 + -) 
u n + ni \ r / 

n + ni V r o / 

The velocity can be directly deduced from first moment of (7) 
11 1 Z" 00 

V = D- <T ™ >= D-J ^r[P 0{ r) +Pl (r)] (23) 



There exists therefore a dynamical phase transition at v = (v\ — i/ ) = 1 

if v < 1 



[r (^i - 1) + ri^o] H ~~ — [to^i + ti(i/ + 1)] 



if > 1 



n + ni n + ni 



(24) 

More generally, the algebraic decay of the distributions Pq,i(t) in the limit r — > oo 

P ,i( r ) oc with ^ = ^-^0 = 1(^1-^°) (25) 

will generate the same succession of anomalous diffusion behaviours for the cumulants 
of the Brownian particle position as a function of the parameter u, as the one that 
appears in the white noise case as a function of the parameter \x [3]. In particular, the 
thermal average of the position of the particle will grow linearly in time only in the 
finite velocity phase 

/+oo 
dx x P(x,t \x Q ,0) oc Vt if v > 1 (26) 

However, in the vanishing velocity phase, the thermal average of the position will grow 
slower than linearly, as a power of time with the exponant v which depends continuously 
on the parameters of the disorder 

/+oo 
dx x P(x,t \x ,Q) oc t v if < v < 1 (27) 

Let us now consider the stationary flux going through a sample between two 

concentrations Po > and Pn = in the limit iV — > oo. The change of variable (10) 

gives immediately the two joint laws Vq(J) et V\{J) from result (21), with the notations 
j = AlQi and j = Dofo 



TO 



V (J)=A9(J) 0(J o -J) J" 1 
Vi(J) = B 6(J) 6(J - J) J" 1 -"" 



, -i (Jq - JY ' 1 



(Ji + jy- 
-i (j -jy° 



(28) 



(Ji + J)^ +1 
where A and B are two normalization constants. 

These two distributions for the flux are concentrated on the bounded interval [0, Jo]. 
The transition at v = v\ — vq = 1 for the self-averaging velocity (24) corresponds for 
the flux probability distributions to a transition of the behaviour in the limit J — > 



V ,i(J) J"' 1 — { } (29) 
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The vanishing velocity phase v = v x — v < 1 is characterised by the divergence of Vq 
and V\ at the origin J — > 0. On the contrary in the finite velocity phase v > 1, the 
probability distributions of the flux vanish at the origin. 

Some curves Vq^{J) are drawn on Figures 2 and 3 for different values of the 
parameters. They illustrate in particular the transitions at v = 1 we just mentioned 
(29) and the transion at vq = 1 for the behaviour of Vo(J) as J — > Jo- We refer to the 
figures in [8] for comparison with the uncorrelated case. 

In this article, we have analysed some properties of anomalous diffusion in a ran- 
dom medium described by a quenched random force {F} distributed as a two-level 
Poisson process. In particular, we showed that there exists a dimensionless parameter 
z/, which governs the asymptotic behaviours of the probability distributions Po,i( r ) m 
the limit r — > oo. As a result, this parameter v also controlls the dynamical phase 
transition for the velocity and the behaviours of the probability distributions "Po,i(J) 
of the flux in the limit J — > 0. All this is qualitatively the same as what is known in 
the white noise case. However, the probability distributions for the random variable 
and for the stationary flux Joo are very different from the white noise case outside the 
asymptotic regimes r — > oo and J — > 0. Indeed, they present restricted supports, as a 
consequence of the bounded character of the random force, and they are given in terms 
of only rational functions and no exponential. 

More generally, the two-level Poisson process we considered is a technically very 
convenient disorder model, since it allows analytical studies despite the presence of 
correlations. It has already been used in various contexts, like one-dimensional quantum 
localization [10] [11], non-linear systems coupled to a random environment [12], and 
noise-induced perturbations on Josephson junctions [13]. The limitation to two-level 
processes simplifies computations, but the approach can be generalized to any finite 
number of levels for the random process [14] . 
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Figures captions 

1) Example of the random potential U(x) seen by the Brownian particle when the 
random force {F} is a two-level Poisson process. 

2) Examples of flux probability distributions Vo(J) (plain line) and "Pi (J) (dashed 
line) in the vanishing velocity phase < v < 1 for 

2a) v o < 1 
2b) u > 1 

3) Examples of flux probability distributions Vo(J) (plain line) and Vi(J) (dashed 
line) in the finite velocity phase v > 1 for 

3a) v o < 1 
3b) v Q > 1 
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